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Inelastic pp collisions are dominated by soft (low momentum transfer) physics where perturba- 
tive QCD cannot be fully applied. A deep understanding of both soft and semi-hard processes is 
crucial for predictions of minimum bias and underlying events of the now coming on line pp Large 
Hadron Collider (LHC). Moreover, the interaction of cosmic ray particles entering in the atmo- 
sphere is extremely sensitive to these soft processes and consequently cannot be formulated from 
first principles. Because of this, air shower analyses strongly rely on hadronic interaction models, 
which extrapolate collider data several orders of magnitude. A comparative study of Monte Carlo 
simulations of pp collisions (at the LHC center-of-mass energy ~ 14 TeV) using the most popular 
hadronic interaction models for ultrahigh energy cosmic ray (SIBYLL and QGSJET) and for col- 
lider physics (the PYTHIA multiparton model) is presented. The most relevant distributions are 
studied including those observables from diffractive events with the aim of discriminating between 
the different models. 



PACS numbers: 13.85.-Tp, 96.40.-z 



I. INTRODUCTION 

Particle colliders and ultra high energy cosmic ray 
(UHECR) experiments are at present the best scenarios 
to test the deep structure of matter and the interactions 
of its fundamental constituents at the frontier of energy. 
To this end, detailed modeling of the underlying physics, 
based on simulation programs also known as event gen- 
erators, are required. These models are very important 
tools to define experimental and analysis strategies, to 
test new theoretical ideas and to design new experiments. 

From the perspective of particle physics, UHECR in- 
teractions are orders of magnitude beyond what can be 
achieved in current (and future) terrestrial collider ex- 
periments and may open a window to energy and kine- 
matic regions previously unexplored in the study of fun- 
damental interactions. The Pierre Auger Collaboration, 
exploding the potential of the hybrid design of the ob- 
servatory (ground based and fluorescence detectors) has 
developed a method to obtain the energy spectrum which 
does not rely on detailed numerical simulations extrapo- 
lated from experimental knowledge of man-made acceler- 
ators to the highest cosmic ray energies [lj]. The hybrid 
detection provides a way to inter-calibrate the subsys- 
tems and to control systematic uncertainties @- This 
new approach to derived the cosmic ray spectrum will 
allow to constrain, although in an indirect manner, inter- 
action models at energies and phase space regions that 
complement those of colliders experiments. For the case 
of primary composition determination, however, UHECR 
data are interpreted in all cases using Monte Carlo sim- 
ulations [111 IE Hi- 

The Large Hadron Collider (LHC) Q , currently under 
construction at CERN, will provide pp collisions at the 
unprecedented centre-of-mass energy of ^fs = 14 TeV 
and luminosity of C = 10 34 cm _2 s _1 . It will also pro- 
vide heavy ion collisions at a centre-of-mass energy of 
about 1000 TeV. A good description of both hard and 
soft processes in event generators at LHC energies are 



very important to understand the underlying event as- 
sociated with high transverse momentum (j?t) jets or 
leptons as well as minimum bias events. These events 
will be responsible for most of the radiation background 
expected at LHC experiments and hence it is essential 
to study detector damages, triggering systems, detector 
counting rates, etc. Minimum bias events are worthy os 
scientific study as they provide a good insight into the in- 
ternal structure of protons. In addition, dedicated runs 
of the LHC with lower luminosity (£ = 10 28 cm _2 s _1 ) 
and specially tuned beam optics are planned to study 
diffractive events. Both ATLAS [1] and CMS[| experi- 
ments are planning to implement additional detectors to 
cover the forward diffractive regions with tracking and/or 
calorimetry [1(J, [III, UM ■ Recently, a new experiment to 
study very forward particle emission in the LHC collider, 
LHCf, has ben approved. The experimental results of 
LHCf will be able to provide the production spectum of 
secondary particles in the very forward region, allowing 
to constrain the Monte Carlo codes Moreover, novel 
aspects of diffraction studies are included in the physics 
case for forward proton tagging at 420 m at the LHC [lj| . 
A study of diffraction must use detectors with excellent 
forward acceptance to allow for a comparison with cos- 
mic ray data. A g ood review of diffractive physics can 
be found in ref . [l5| |. while references [IE [IE EE OH focus 
on future studies at LHC. 

In this paper we present a comparative study between 
the two most frequently applied models for simulation 
of cosmic rays extensive air showers SIBYLL [20j and 
QGSJET 12111 with a multi-purpose Monte Carlo like 
PYTHIA [22| tunned for use in LHC experiments. The 
paper is organized as follows: In section II the main fea- 
tures of the models used in the Monte Carlo generators 
for pp collisions with emphasis in the most distinctive 
differences among them are presented. The method and 
the results from the Monte Carlo studies at LHC center 
of mass energy with the aim of discriminating between 
the different models are described in section III with two 
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parts. In the first part, predictions for the most relevant 
distributions in studies of both collider and UHECR data 
are discussed. In the second part a comparative analy- 
sis of the signatures from diffractive events is presented. 
Final comments and conclusions are given in section IV. 



II. MODELS FOR HADRONIC COLLISIONS 

Although electromagnetic and weak interactions are 
well understood, this is not the case for hadron pro- 
duction in collisions of nucleons, pions and kaons with 
light nuclei, where the lack of experimental data posses 
limitations in many UHECR and accelerator applica- 
tions [J, |23| . This is mainly because precise pertur- 
bative quantum chromo-dynamics (pQCD) calculations 
are only possible for processes with large momentum 
transfer, also known as "hard" interactions, that consti- 
tute only a minute fraction of the overall reaction rate. 
In hadron-hadron interactions it is customary to distin- 
guish between elastic and inelastic processes, and this 
later into diffractive (including single and double diffrac- 
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tion) and non-diffractive ones (usually called minimum 
bias events). Precisely, low-px ("soft") processes, where 
pQCD cannot be fully applied and phenomenology mod- 
els are used, play a dominant role in the non-diffractive 
component. 

Current models of high energy hadron collisions typi- 
cally rely on the pQCD formalism in the description of 
high-j»T scattering, while treating the low-pT ones in a 
simplified phenomenological approach. At the LHC, the 
description of the new physics processes to be studied are 
mostly controlled by pQCD. Some efforts have been done 
to investigate the models used by the most popular col- 
lider event generators, like PYTHIA, to describe "soft" 
interactions in hadron-hadron collisions with the aim of 
predicting minimum bias and underlyin g ev ent levels of 
particle production at the LHC [H [1| H|, H3, H HH • 
These studies resulted in a tuning of PYTHIA based on 
comparisons of experimental data that is used in this pa- 
per. 

The pQCD inclusive cross section of production of par- 
ton jets pairs with transverse momenta larger than some 
cutoff Qf nin is given by, 



vqcb(s,pT° S )=J2 
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where x\ and x 2 are the fractions of the momenta of 
the parent hadrons carried by the partons which col- 
lide, d<Tij/d\i\ is the cross section for scattering of par- 
tons of types i and j according to elementary QCD 
diagrams, fi and fj are parton distribution functions 
(pdf's), § = x\ X2S and — t = s(l — cos$*)/2 = Q 2 are 
the Mandelstam variables for this parton-parton process, 
and the sum is over all parton species. 

In the UHECR field, the required information to 
model the interaction of the primary particle entering 
in the atmosphere appears to be extremely sensitive to 
the underlying "soft" non perturbative hadronic process 
[3p| . In this direction, there are three event generators, 
SIBYLL [H, QGSJETdH and DPMJET[3l| which are 
tailored specifically for simulation of hadronic interac- 
tions up to the highest cosmic ray energies. 

The most frequently used high energy hadronic mod- 
els in the study of UHECRs are QGSJET and SIBYLL. 
In these codes, the low pr interactions are modeled by 
the exchange of Pomerons. Regge singularities are used 
to determine the momentum distribution functions of 
the various sets of constituents, valence and sea quarks. 
Both QGSJET and SIBYLL share the eikonal model and 
then assume the unitarized cross sections assuming a real 



eikonal function sum of a soft and hard contributions: 



OmeZ = / d 2 b 



*(l - cxp{-2 Xs (s, b) - 2 Xh (s, &)}) (2) 



At high energies the hard eikonal is dominating: 



(3) 



where the normalized profile function A(s, b) describes 
the distribution of partons in the plane transverse of the 
collision axis. QGSJET and SIBYLL take different as- 
sumptions on the profile function which determines the 
inelastic cross section and its energy dependence. 
QGSJET assumes a gaussian profile distribution and its 
theory is formulated entirely in terms of Pomeron ex- 
changes. The basic idea is to replace the soft Pomeron 
by a so-called semihard Pomeron, defined to be an ordi- 
nary soft Pomeron with the middle piece replaced by a 
QCD parton ladder. Thus, minijets will emerge as a part 
of the semihard Pomeron, which is itself the controlling 
mechanism for the whole interaction. 

In SIBYLL the profile function is based on the Fourier 
transform of the electromagnetic form factor and it is 
an energy-independent exponential. The underlying idea 
behind SIBYLL is that the increase in the cross section 
is driven by the production of minijets. The probability 
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distribution for obtaining N jet pairs (with pj, > p™ ln , 
being p™ ln a sharp threshold on the transverse momen- 
tum below which hard interactions are neglected) in a 
collision at energy y/s is computed regarding elastic pp or 
pp scattering as a diffractive shadow scattering associated 
with inelastic processes. The algorithms are tuned to re- 
produce the central and fragmentation regions data up 
to pp collider energies, and with no further adjustments 
they are extrapolated several orders of magnitude. 

A general update of QGSJET has been recently pre- 
sented, where the key improvement is connected to an 
account for non-linear interaction effects in individual 
hadronic collisions [HJ . Additionally, a more reliable low 
mass diffraction treatment has been used and all model 
parameters have been re-calibrated using a wider set of 
accelerator data. This results in a new model, QGSJET- 
II. Non-linear screening corrections appear to be corre- 
lated with corresponding parton densities and become 
larger at higher energies, smaller impact parameters, re- 
sulting in the saturation of pdf's at the scale Q^ ain and 
in a considerable reduction of "soft" particle production. 

In the case of PYTHIA, perturbative QCD is used ex- 
tending it for the case of low-pt. pQCD is divergent for 
Pt 0, PYTHIA avoids the divergence using two dif- 
ferent scenarios. The "simple scenario" consists in fixing 
a minimum value of p™ m below which the cross section 
is defined as null and can be interpreted as the inverse of 
some color screening length in the hadron. This is equiv- 
alent to set a maximum impact parameter 6 max above 
which there is no more interaction. In the so called "com- 
plex scenario" a regulating parameter ptq is introduced 
below which the cross sections are dumped. Different 
models of matter distribution in the hadron are consid- 
ered: uniform, simple gaussian and double gaussian. 

The transition process from asymptotically free par- 
tons to colour-neutral hadrons is describes in all codes 
by string fragmentation models [34| . 

In summary, there are differences between the models 
for hadronic collisions in the existing event generators 
that will emerge in the Monte Carlo study presented in 
the rest of this paper. 



III. DESCRIPTION OF THE METHOD AND 
MONTE CARLO STUDIES 

For analyzing the difference between models, we have 
generated samples of 10 4 pp collisions at the LHC 
center-of-mass energy for each event generator model: 
QGSJET-II [H, QGSJET-01 [H, SIBYLL 2.1 [H and 
PYTHIA 6.205 [221] . All calculations contain a mixture 
of diffractive and non-diffractive events according to the 
model used. All secondary particles were registered with- 
out any energy cut. 

In the case of PYTHIA, it was chosen the "complex 
scenario" with a double gaussian distribution of matter 
inside the hadron. The default values of some parame- 
ters were modified according to the results of references 
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FIG. 1: Distribution of the number of charged secondaries 
produced in pp collision at LHC energy. The solid line cor- 
responds to PYTHIA 6.2, dotted to QGSJET-01, the dashed 
one to SIBYLL2.1 and dashed-dotted to QGSJET-II 



0, [Mill, IIS [11, III where the optimal values where ob- 
tained from a tunned PYTHIA using events from differ- 
ent experiments. In the following table both default and 
tunned values of PYTHIA as used in this paper are pre- 
sented: where MSTP(81) refers to the master switch for 



TABLE I: Values of relevant PYTHIA 6.2 parameters. 



Variable 


Default 


Tunned 


MSTP(81) 


1 


1 


MSTP(82) 


1 


4 


PARP(82) 


2.1 


1.8 


PARP(83) 


0.5 


0.5 


PARP(84) 


0.2 


0.5 


MSTP(2) 


1 


1 



multiple interactions, MSTP(82) = 4 selects an hadronic 
matter overlap consistent with a given double Gaussian 
matter distribution and a continuous turn-off of the cross 
section at pto =PARP(82). This double Gaussian mat- 
ter distribution is regulated by the following parameters: 
a core PARP(84) of the main radius containing a frac- 
tion PARP(83) of the total hadronic matter. The value 
of the parameter MSTP(2) gives the kind of calculation 
of as at hard interaction and if the value is 1 then it is 
first-order running of as (Here as is the strong coupling 
constant) 
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TABLE II: Average multiplicity and inelasticity per proton- 
proton collision. 
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80.3 
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76.1 


N total 


126.5 


139.3 
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<k L > 


0.41 


0.50 


0.43 


0.43 
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To help understanding the differences of the models 
currently used in the study of UHECRs, when extrapo- 
lated from collider data to higher energies and to small 
angle processes, it is important to compare those vari- 
ables driving the development of air showers, i.e, the 
multiplicity of produced secondary particles and the in- 
elasticity (the relative energy loss of leading secondaries) . 

Multiplicity distributions of charged secondaries (N^) 
produced in each collision are shown in Fig[TJ Table 2 
shows the average number of secondaries for each model. 
Besides the mean values, the standard deviations over 
the 10 4 interactions are given. The errors of the mean 
values (er/-\/(10 4 )) are much smaller. SIBYLL produces 
60% to 65% less nucleons than the other models. This 
has been noticed in previous analysis at higher energies 
@, H and it is also seen at LHC energy. QGSJET-01 
produces more baryons and charged pions than the other 
models, but this situation has changed in QGSJET-II, 
due to the non- linear screening corrections, which lead 
to a reduction of the interaction eikonal and hence of 
the number of elementary particle production processes 
f32| . This results in a reduction of particle produc- 
tion in QGSJET-II compared with QGSJET-01, bring- 
ing it closer to SIBYLL and in good agreement with the 
PYTHIA predictions. The mean charged particle mul- 
tiplicity, which already shows differences between mod- 
els at this energy, readily increases with rising energy as 
QGSJET predicts a power law-like increase of the num- 
ber of secondaries, while SIBYLL multiplicity exhibits a 
logarithmic growth and PYTHIA follows a ln 2 (s). It is 
worth mentioning here the results reported by the CDF 
collaboration favouring an energy dependence stronger 
than ln(s) [H- 

The distribution in pseudorapidity, r\ = — lntan(#/2), 
of charged particles for pp collisions at 14 TeV is pre- 
sented in Fig. H It shows up clearly that QGSJET-01 
produces more secondaries than the other models in all 
directions but in the central region of |r?|, where it is 
superseded by PYTHIA. This might be explained by the 



FIG. 2: Pseudorapidity distribution for charged particles. 



pure QCD treatment and the possibility of multiple inter- 
actions set in PYTHIA. It has been noticed [25[ that this 
tunned version of PYTHIA provides the best description 
of experimental data from UA5 and CDF in the central 
rapidity region. Above \r)\ > 3 PYTHIA and SIBYLL 
give similar predictions, both having smaller values than 
the QGSJET models. 

To get additional information, the pseudorapidity dis- 
tributions for charged and neutral pions, kaons, protons, 
antiprotons and neutrons are shown in Fig. [3J There 
are discrepancies between the models in all cases being 
the largest ones for nucleons and antiprotons, where al- 
though models predict similar shapes, SIBYLL exhibits 
a clear deficit at all rj. For the case of kaons there is a 
factor of two in the predictions from the two versions of 
QGSJET. The peaks in the very forward and backward 
parts of the pseudorapidity distribution for protons corre- 
spond to diffractive events, in which one of the smashing 
protons keeps traveling approximately in the same direc- 
tion after the collision. The shoulders in the high pseu- 
dorapidy region for neutrons have a different origin. If a 
neutron comes out as the fast particle, charge is being ex- 
changed. The process, accordingly, cannot be attributed 
to the exchange of zero quantum numbers (i.e to diffrac- 
tion) but, for instance, to pion exchange. It is worth 
noting that the differences in the production of neutral 
pions influence the shower development of the secondary 
particles produced by the interaction of a primary cos- 
mic ray particle in the atmosphere, which is driven by 
the electromagnetic component generated from the 7r°s. 

Two-dimension distributions of number of secondaries, 
N sec , vs pseudorapidity, i], are presented in Fig |U 
The plots show in detail the large differences between 
QGSJET, upper panels, and SIBYLL, bottom-left panel, 
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FIG. 3: Pseudorapidity distribution for different types of par- 
ticles. Left-top panel corresponds to charged pions, right- 
top to neutral pions, left-bottom to protons and right-bottom 
panel corresponds to neutrons. The solid line corresponds to 
PYTHIA 6.2 events, dotted to QGSJET-01, the dashed one 
to SIBYLL2.1 events, dashed-dotted to QGSJET-II. 



in the whole region of rj. The two diffractive peaks in 
the region of low multiplicity and high pseudorapidity 
are well separated from the broad distribution of non- 
diffractive events in both versions of QGSJET and to 
less extent in PYTHIA. This feature is due to the fact 
that QGSJET models have none or few non-diffractive 
events with small multiplicities, while in PYTHIA the 
distributions overlap and diffractive events tend to have 
higher multiplicities as well. For the SIBYLL model the 
distribution for low number of secondaries is rather flat 
in all directions. 

Table III shows the percentage frequency of the leading 
particle produced in the collision. In more than 50% of 
the collisions protons emerge as these leading particles. 
In SIBYLL an almost 65 % the most energetic particle 
is a proton, while the other models give between 43% 
to 62%. SIBYLL and PYTHIA generate mesons as the 
most energetic particle in ss 20% of the cases, however 
QGSJET-01 and QGSJETII have larger and smaller pro- 
ductions of mesons as leading particles respectively with 
a difference of roughly a factor of 3. All models assume 
that the leading particle distributions scale with energy, 
being tunned to low energy. Certainly, measurements of 
hadron production in the very forward region are needed 
to study the leading baryon distributions, mainly because 
there are some theoretical models predicting that the 
leading particle distributions will change drastically at 
very high energies [H, H3] 



FIG. 4: Number of secondaries vs pseudorapidity 2D- 
distribution. Top-left panelcorresponds to QGSJET-II, top- 
right panel corresponds to QGSJET-01, bottom-left panel to 
SIBYLL 2.1 and bottom-right panel to PYTHIA 6.2. 



TABLE III: Most energetic secondary particle probabilities. 





Pythia 6.2 


Qgsjet 01 


Qgsjet II 


Sibyll 2.1 


proton 


55.29% 


43.27% 


62.08% 


64.62% 


neutron 


27.34% 


18.31% 


19.68% 


16.51% 


£ nucleons 


82.63% 


61.58% 


78.76% 


81.13% 


7T± 


10.28% 


20.47% 


6.77% 


10.72% 


7T° 


4.89% 


9.74% 


3.02% 


5.85% 




1.57% 


2.40% 


0.73% 


1.00% 


K L 


0.63% 


0.91% 


0.44% 


0.57% 



B. Signatures of diffractive events 

As mentioned above, in hadron-hadron interactions the 
inelastic processes are usually divided into diffractive and 
non-diffractive. In this section, a study of the predicted 
signatures from different models for diffractive hadronic 
interactions is presented. 

A good parameter for tangling diffractive events from 
pp collisions is the inelasticity defined as: 



where Ep is the energy of the incident particle in the 
lab frame, and -Eiead is the energy of the secondary with 
largest energy (the so-called leading particle). A signa- 
ture that can be used to distinguish diffractive from non- 
diffractive events is their low value in both inelasticity 
and number of secondaries [38[ . In Table II the value of 
the average inelasticity for each model is shown, while 
the corresponding inelasticity distributions are displayed 
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FIG. 5: Inelasticity (fcr = 1 - %^-) distriution. The solid 
line corresponds to PYTHIA 6.2, dotted line to QGSJET- 
01, dashed line to SIBYLL2.1 and dashed-dotted one to 
QGSJET-II . 



FIG. 6: Distribution of the number of secondaries for events 
with low inelasticity (fc^ < 0.04). 



in Fig. [5J A narrow peak at low fcz, is evident from 
this plot, which corresponds to elastic and single diffrac- 
tive processes. For non-diffractive events, the available 
energy is shared among many secondaries leading to a 
rather uniform distribution in the whole range of k^. 

A few differences between models can be seen in the in- 
elasticity distribution: PYTHIA has the highest diffrac- 
tive peak. QGSJET-01 also presents a small peak at large 
ki due to the large number of secondaries produced in 
inelastic collisions. QGSJET-II has no longer that fea- 
ture. SIBYLL and PYTHIA distributions are in general 
good agreement. 

A close inspection of the multiplicity for single diffrac- 
tive, double diffractive and non-diffractive events using 
PYTHIA indicates that cutting at multiplicity below 
40 the diffractive events largely dominates the sample. 
However, a cut in low inelasticity (for, < 0.04) leave 
an even more pure diffractive sample. We will then 
label as "diffractive" events the ones with inelasticity 
< 0.04. In Fig. [5] the distribution of charged par- 
ticles, N ch , for k L < 0.04 is presented. Both QGSJET- 
01 and SIBYLL distributions barely goes further than 
40 secondaries, while QGSJET-II and PYTHIA extend 
up to 60 secondaries. Figure [7] shows particle densi- 
ties distributed in pseudorapidity space for "diffractive" 
events. There are large divergencies between PYTHIA 
and SIBYLL in the predictions of the particle multiplic- 
ity in the whole region of rj. PYTHIA predicts a density 
roughly 80 % greater than SIBYLL. QGSJET-01 shows 
a rather flat distribution at intermediate values between 
PYTHIA and SIBYLL, while the new version QGSJET- 
II is in good agreement with PYTHIA at < 4 and 



present smaller values at \r)\ > 4. It is worth mention- 
ing here that the cuts kh < 0.04 and N sec < 40 in the 
PYTHIA sample allow selection of a pure sample of single 
diffractive events. 

The transverse momentum, px, detection capabilities 
at ATLAS and CMS with good resolution will be lim- 
ited to particles with p T > 0.5 GeV @, @ where for the 
labeled "diffractive" events QGSJET-II and PYTHIA 
show the largest difference. This is evident in Fig. [8] 
where the pt distribution of eventsin the central region ( 
\rj\ < 5) is shown. At low momenta dN^/dpr is greater 
for PYTHIA, but as pr increases densities for the other 
models become greater with a difference of an order of 
magnitud at px > 1.5 for the case of QGSJET-II, as 
PYTHIA does not create diffractive events with high 
transverse momemtum. 

Figure [5] shows 2-D distributions of the number of sec- 
ondaries in the labeled "diffractive" events vs. pseudora- 
pidity. Again, it is evident that the distributions gener- 
ated by SIBYLL are fundamentally different to the ones 
from QGSJETs and PYTHIA. This is certainly due to 
the phenomenological description of diffractive events in 
SIBYLL [H. QGSJET-01 and QGSJET - II predict a 
large fraction of events with low number of secondaries 
in the region of high pseudorapidity. Clearly, the better 
treatment of diffraction for the case of QGSJET-II with 
its parameters tunned using accelerator data, results in 
diffractive peaks much lower than in QGSJET-01, but 
still larger compared with PYTHIA data. 
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FIG. 7: Pseudorapidity distribution for the selected "diffrac- 
tive" events. 
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FIG. 8: pr distribution for the selected "diffractive" events. 



IV. CONCLUSIONS 

The two most widely used high energy hadronic inter- 
action models for the study of cosmic rays, SIBYLL and 
QGSJET, and the most popular simulation program for 
collider experiments, PYTHIA, have been compared at 
LHC center of mass energy. These codes using different 
theoretical models involving many variables have differ- 
ent predictions for the most important observables, such 
as particle multiplicities, distribution of particles in pseu- 



FIG. 9: Multiplicity vs pseudorapidity 2D-distribution for low 
inelasticity (fez, < 0.04) events. Top-left panel corresponds 
to QGSJET-II, top-right panel corresponds to QGSJET-0I , 
bottom- left panel to SIBYLL 2.1 and bottom-right panel to 
PYTHIA 6.2. 



dorapidity space and transverse momentum distribution, 
allowing to investigate, and improve, the theoretical pre- 
dictions of hadron-hadron interactions at this energy. 

The discrepancies observed in the pseudorapidity dis- 
tributions for different types of particle, in particular in 
the acceptance region of LHC experiments would allow 
to perform a combined analysis to constrain the models, 
using the large minimum bias events statistics that will 
be collected in the very first LHC operation. The predic- 
tions for the 2-D distributions of particle multiplicity vs 
pseudorapidity will also help deciding on the best theo- 
retical model. 

Measurements of hadron production with the forward 
detectors attached to LHC experiments are very impor- 
tant for a better understanding of the simulations that 
model soft hadronic interactions at high energies. The 
analysis of the percentage frequency of the leading par- 
ticle produced in the collision indicates differences up to 
a factor of 2 for meson production. 

A selection of events from pp collisions at *Js = 14TeV 
with small inelasticity (k^ < 0.04) and low number of 
secondaries allows to pick diffractive events for a compar- 
ative study of the various models. The analysis of particle 
densities in pseudorapidity space indicates a good agree- 
ment of predictions using PYTHIA and QGSJET-II and 
a clear deficit of particle densities using SIBYLL and to 
less extent QGSJET-01. The pr distributions in the cen- 
tral pseudorapidity region clearly indicate that SYBYLL 
and QGSJET models create up to an order of magnitude 
more particles with large pt than PYTHIA. A compari- 
son with experimental data will provide strong constrains 
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in modeling diffractive physics. 

The lack of suitable accelerator data is the dominant 
source of systematic uncertainties in the analysis of the 
extensive air shower data. At the same time, analysis 
of minimum bias events at LHC are very important for 
understanding the underlying event and commissioning 
studies for LHC detectors. Certainly, the discrepancies 
in the models discussed in this paper will be naturally 
reduce with the large statistics of interesting data at a 
completely new energy frontier for terrestrial colliders, as 
LHC, and cosmic ray experiments, as the Pierre Auger 
Observatory. 
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